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RESERVOIR MODEL BUILDING METHODS 



CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] application claims the benefit of U.S. Provisional Application 

60/540,794, filed January 30, 2004. 

BACKGROUND 
Field of Invention 

[0002] Embodiments of the invention are related to evaluation of subsurface 
reservoirs. 

Description of Related Art 

[0003] ^ industry, geologic models are often used to aid in 

activities, such as determining the locations of wells, estimating hydrocarbon 
reserves, or plaiming reservoir-development strategies, including evaluating the 
economic recovery of hydrocarbon resources. A geologic model typically is a 
computer-based representation of a subsurface earth volume, such as a petroleum 
reservoir or a depositional basin. 

[0004] Geologic models may take on many different forms. Depending on the 
context, descriptive or static geologic models built for petroleum applications can be 
in the form of a 3-D array of cells, to which geologic and/or geophysical properties 
such as lithology, porosity, acoustic impedance, permeability, or water saturation are 
assigned (such properties will be referred to collectively herein as "reservoir 
properties"). 

[0005] Many geologic models are constrained by stratigraphic or structural 
siurfaces (e.g., flooding surfaces, sequence interfaces, fluid contacts, faults) and 
boundaries (e.g., facies changes). These surfaces and boimdaries define regions 
within the model that possibly have different reservoir properties. 
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[0006] Various approaches can be followed for evaluating a reservoir using 
geologic modeling. At least one approach is strictly sequential, involving sequential 
evaluations by several disciplines. With such an approach, a reservoir evaluation 
using geologic modeling might take several or many months to complete. With such 
an approach, due to the large amount of time necessary for evaluating a reservoir 
using geologic modeling, only one geologic model would tend to be built in 
connection with the reservoir evaluation. Consequently, such an approach would 
allow no realistic opportunity to leam how decisions are made during the geologic 
modeling process, or how such decisions would affect the final outcome. Such a strict 
sequential approach would also allow no opportunity to evaluate the inherent 
imcertainty in arriving at solutions to problems, considering the limited amoimt of 
data that would tend to be available for use in the geologic modeling as well as the 
level of interpretation required in the geologic modeling process. 

[0007] Furttiermore, such a strict sequential £^proach for evaluating a reservoir 
using geologic modeling would in all likelihood tend to involve building a geologic 
model made up of many millions of cells, e.g., 200 million cells, and require 
**upscaling'* the geologic model in order to reduce the niunber of cells to no more than 
500,000 cells so that flow simulation could be performed. Obviously, the steps of 
building geologic models and then upscaling them would tend to contribute further to 
the large amoimt of time needed to evaluate a reservoir using geologic modeling. 

[0008] Accordingly, a need exists for improved methods of evaluating reservoir. 
SUMMARY 

[0009] Embodiments of the invention are directed to a method for generating a 
reservoir model. In one embodiment, the method includes providing a first 
firamework having a plurality of cells, wherein the first fi-amework is a reservoir 
fi-amework and providing a second fi-amework having a plurality of cells, wherein the 
volume of the first fi-amework is greater than the volume of the second firamework. 

[0010] another embodiment, the method includes providing a first firamework 

having a plurality of cells, wherein the first firamework is a reservoir firamework, and 
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providing a second framework having a plurality of cells, wherein the volume of the 
second framework is substantially the same size as one of the cells of the first 
framework. 

[001 1] In yet another embodiment, the method includes providing a fi-amework 
having a plurality of cells, wherein each cell is substantially the same size as a sample 
of the well data, identifying some or all of the cells of the framework as net or non- 
net, populating some or all of the cells of the framework with one or more reservoir 
properties to provide a reservoir cell model, and performing a flow simulation on the 
reservoir cell model to generate one or more effective reservoir property values. 

[0012] ^ still another embodiment, the method includes providing a first 
framework having a plurality of cells, wherein the first framework is a reservoir 
framework, and providing a second framework having a plurality of cells, wherein 
each one of the cells of the second framework is substantially the same size as a 
sample of the well data. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0013] Figure 1 illustrates a flow diagram of a method for generating one or more 
effective reservoir properties in accordance with one embodiment of the invention. 

[0014] Figure 2 illustrates a cell framework in accordance with one embodiment 
of the invention. 

[0015] Figure 3 illustrates a flow diagram of a method for generating one or more 
reservoir property values in accordance with another embodiment of the invention. 

[0016] Figure 4 illustrates a computer network, into which embodiments of the 
invention may be implemented. 

DETAILED DESCRIPTION 
Introduction and Definitions 

[0017] ^ detailed description will now be provided. Each of the appended claims 
defines a separate invention, which for infringement purposes is recognized as 
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including equivalents to the various elements or limitations specified in the claims. 
Depending on the context, all references below to the "invention" may in some cases 
refer to certain specific embodiments only. In other cases it will be recognized that 
references to the "invention" will refer to subject matter recited in one or more, but 
not necessarily all, of the claims. Each of the inventions will now be described in 
greater detail below, including specific embodiments, versions and examples, but the 
inventions are not limited to these embodiments, versions or examples, which are 
included to enable a person having ordinary skill in the art to make and use the 
inventions, when the information in this patent is combined with available 
information and technology. Various terms as used herein are defined below. To the 
extent a term used in a claim is not defined below, it should be given the broadest 
definition persons in the pertinent art have given that term as reflected in printed 
publications and issued patents. 

[0018] ^ herein, the term "cell" is defined as a unit or block that defines a 
portion of a three dimensional reservoir model. As such, a three dimensional 
reservoir model may include a number of cells, ranging Scorn tens and hundreds to 
thousands and millions of cells. Each cell represents a specifically allocated portion 
of the three dimensional reservoir model. An entire set of cells may constitute a 
geologic model and thus represent the subsurface earth volume of interest. Each cell 
preferably represents a unique portion of the subsurface. As such, the cells preferably 
do not overlap each other. Dimensions of the cells are preferably chosen so that the 
reservoir properties within a cell are relatively homogeneous, yet without creating an 
excessive number of cells. Preferably, each cell is square or rectangular in plan view 
and has a thickness that is either constant or variable. However, it is contemplated that 
other shapes may alternatively be used. 

[0019] ^ used herein, the term "reservoir properties" are defined as quantities 
representing physical attributes of rocks containing reservoir fluids. The term 
"reservoir properties" as used in this application includes both measurable and 
descriptive attributes. Examples of measurable reservoir property values include 
rock-type firaction (e.g., net-to-gross, v-shale, or facies proportion), porosity, 
permeability, water saturation, and fi^cture density. Examples of descriptive reservoir 
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property values include facies, lithology (e.g. sandstone or carbonate), and 
environment-of-deposition (EOD). Reservoir properties may be populated into a 
reservoir framework to generate a reservoir model. 

[0020] The term "rock-type fraction" is defined as the ratio of the rock volvmie 
containing a specific rock-type that to the total (gross) rock volxmie. As such, the 
gross rock volume can be divided into 2 components: (1) rock volume containing a 
specific rock-type, and (2) rock volume containing all other rock types. So, rock-type 
fi-action may be expressed as: 

volume of a specific rock-type 

rock-type firaction = 

total rock volume 

[0021] Example of a rock-type fraction is v-shale (volimie shale), typically 
calculated from electronic well log measurements and sometimes inferred from 
seismic data. Using the expression for rock-type fraction: 

volume of shale 

v-shale = 

total rock volume 

[0022] The term "net-to-gross", also denoted N:G, as used herein includes the 
term v-shale (volume shale). The relationship between v-shale and net-to-gross may 
be expressed as follows: 

net-to-gross = 1 - v-shale. 

Furthermore, whenever the term "net-to-gross" or **N:G" is used herein, it will be 
understood that this is an example of a rock-type fiaction, and that any other choice of 
rock-type fi^ction may be selected. 

[0023] As used herein, the term "permeability" is defined as the ability of a rock 
to transmit fluids through interconnected pores in the rock. Permeability can vary 
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substantially within a hydrocarbon-bearing reservoir. Typically, permeabilities are 
generated for fine-scale models (geologic models) using data from well core samples. 
For simulation cells, the heterogeneities of the geologic model are accounted for by 
determining an effective permeability. An effective permeability of a heterogeneous 
medium is defined as the permeability of an equivalent homogeneous medium that, 
for the same boundary conditions, would give the same flux (amount of fluid flow 
across a given area per unit time). 

[0024] As used herein, the term "porosity" is defined as the percent volume of 
pore space in a rock. Porosity is a measure of the reservoir rock's storage capacity for 
fluids. Porosity is preferably determined firom cores, sonic logs, density logs, neutron 
logs or resistivity logs. Total or absolute porosity includes all the pore spaces, 
whereas effective porosity includes only the interconnected pores. 

[0025] As used herein, the term **well data" is defined as any data that may be 
obtained firom a well. Well data include, but are not limited to, log data and core data. 

[0026] As used herein, the term "geostatistical estimation" is defined as a 
statistical estimon technique used to spatially correlate random variables in geological 
or geophysical applications. Geostatistical estimation involves techniques for 
interpolation and extrapolation of physical measurements using correlation and 
probability concepts. More specifically, geostatistical estimation takes into account 
distance, direction, and spatial continuity of the reservoir property being modeled. 
Geostatistical estimation may be either deterministic or probabilistic. Deterministic 
geostatistical estimation calculates a minimum-variance estimate of the reservoir 
property at each cell. Probabilistic geostatistical estimation develops distributions of 
the reservoir property values and produces a suite of geologic models for the reservoir 
property being modeled, with each model theoretically being equally probable. The 
spatial continuity of a reservoir property may be captured by a variogram, a well- 
known technique for quantifying the variability of a reservoir property as a function 
of separation distance and direction. 
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[0027] As used herein, the term **flow simulation" is defined as a numerical 
method of simulating the transport of mass (typically fluids, such as oil, water and 
gas), energy, and momentum through a physical system using a computer. The 
physical system includes a three dimensional reservoir model, fluid properties, the 
number and locations of wells. Flow simulations also require a strategy (often called 
a well-management strategy) for controlling injection and production rates. These 
strategies are typically used to maintain reservoir pressure by replacing produced 
fluids with injected fluids (e.g. water and/or gas). When a flow simulation correctly 
recreates a past reservoir performance, it is said to be "history matched," and a higher 
degree of confidence is placed in its ability to predict the fixture fluid behavior in the 
reservoir. 

[0028] As used herein, the term *three dimensional reservoir model" is defined as 
a three dimensional firamework of cells that contain reservoir property values. 

[0029] As used herein, the term "three dimensional fi-amework" is defined as a 
numerical representation of a volume that is divided into cells. The numerical 
representation includes the total number of cells, their dimensions, and how they are 
connected to each other. 

[0030] As used herein, the term "target reservoir fi-amework" refers to the 
firamework for a target reservoir model. 

[0031] As used herein, the term "target reservoir model" is defined as a target 
reservoir fi-amework populated with reservoir properties. The target reservoir model 
may be any reservoir model, such as a geologic model, flow simulation model, and 
the like. 

Specific Embodiments 

[0032] Various specific embodiments are described below, at least some of which 
are also recited in the claims. 

[0033] In at least one specific embodiment, a method of generating a reservoir 
model, includes: providing a first framework having a pluraUty of cells, wherein the 
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first fiamework is a reservoir firamework; and providing a second framework having a 
plurality of cells, wherein the volume of the first fi-amework is greater than the 
volimie of the second firamework. 

[0034] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the volume of the second firamework is substantially the 
same size as one of the cells of the first framework. 

[0035] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, each one of the cells of the second firamework is 
substantially the same size as a sample of well data. 

[0036] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, each one of the cells of the second framework is 
substantially the same size as a sample of core data. 

[0037] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, each one of the cells of the second framework is 
substantially the same size as a sample of log data. 

[0038] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiurther includes identifying some or all of the cells of the 
second framework as net or non-net. 

[0039] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes identifying some or all of the cells of the 
second framework as sand or shale. 

[0040] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes fiirther includes populating some or all of 
the cells of the second framework with net and non-net values. 

[0041] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes receiving one or more estimated net-to- 
gross values of the first framework. 
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[0042] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes receiving one or more estimated net-to- 
gross values of the first fi-amework; and identifying some or all of the cells of the 
second firamework as net or non-net according to the estimated net-to-gross values of 
the first framework, 

[0043] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes receiving one or more estimated net-to- 
gross values of the first framework; and populating some or all of the cells of the 
second framework with net and non-net values according to the estimated net-to-gross 
values of the first firamework. 

[0044] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes populating some or all of the cells of the 
second framework with one or more reservoir property values. 

[0045] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes populating some or all of the cells of the 
second framework with one or more porosity values. 

[0046] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes populating some or all of the cells of the 
second framework with one or more permeability values. 

[0047] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes populating some or all of tiie cells of the 
second framework with one or more water saturation values. 

[0048] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes populating some or all of the cells of the 
second framework with one or more reservoir property values to generate a reservoir 
cell model; and performing a flow simulation on the reservoir cell model to generate 
one or more effective reservoir property values for the first framework. 
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[0049] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: populating some or all of the cells of the 
second fi-amework with one or more reservoir property values to generate a reservoir 
cell model; performing a flow simulation on the reservoir cell model to generate one 
or more effective reservoir property values for the first fi-amework; and calculating the 
variability between the effective reservoir property values for the first firamework. 

[0050] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiulher includes: populating some or all of the cells of the 
second framework with one or more reservoir property values to generate a reservoir 
cell model; performing a flow simulation on the reservoir cell model to generate one 
or more effective reservoir property values for the first firework; calculating the 
variability between the effective reservoir property values for the first firamework; and 
determining whether the rate of change in the variability between the effective 
reservoir property values remains substantially the same. 

[0051] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: populating some or all of the cells of the 
second firamework with one or more reservoir property values to generate a reservoir 
cell model; performing a flow simulation on the reservoir cell model to generate one 
or more effective reservoir property values for the first firamework; and populating the 
first firework with the effective reservoir property values to generate the reservoir 
model. 

[0052] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the reservoir model is a flow simulation model. 

[0053] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the reservoir model is a geologic model. 

[0054] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the volume of the second framework is greater than the 
size of one cell of the first framework. 
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[0055] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the second fi-amework includes two or more cell samples 
of the first firamework, wherein each cell sample is substantially the same size as one 
of the cells of the first firamework. 

[0056] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: populating some or all of the cells of the 
second firamework with one or more reservoir property values to generate a reservoir 
cell model; and extracting one or more cell samples firom the reservoir cell model, 
wherein each cell sample is substantially the same size as one of the cells of the first 
fi-amework. 

[0057J A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiuther includes: populating some or all of the cells of the 
second firamework with one or more reservoir property values to generate a reservoir 
cell model; and extracting one or more cell samples firom the reservoir cell model, 
wherein each cell sample is substantially the same size as one cell of the first 
fi-amework; and performing a flow simulation on the cell sample to generate one or 
more effective reservoir property values. 

[0058] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the second framework is three dimensional. 

[0059] In at least one specific embodiment, a method of generating a reservoir 
model, includes: providing a first framework having a plurality of cells, wherein the 
first framework is a reservoir framework; and providing a second framework having a 
plurality of cells, wherein the volume of the second framework is substantially the 
same size as one of the cells of the first framework. 



[0060] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, each one of the cells of the second fi-amework is 
substantially the same size as a sample of well data. 
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[0061] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, each one of the ceils of the second firamework is 
substantially the same size as a sample of core data. 

[0062] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, each one of the cells of the second firamework is 
substantially the same size as a sample of log data. 

[0063] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: identifying some or all of the cells of the 
second fi'amework as net or non-net, 

[0064] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: identifying some or all of the cells of the 
second firamework as sand or shale. 

[0065] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes: populating some or all of the cells of the 
second firamework with net and non-net values. 

[0066] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: receiving one or more estimated net-to- 
gross values of the first firamework. 

[0067] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes: receiving one or more estimated net-to- 
gross values of the first framework; and identifying some or all of the cells of the 
second framework as net or non-net according to the estimated net-to-gross values of 
the first firamework. 

[0068] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes: receiving one or more estimated net-to- 
gross values of the first framework; and populating some or all of the cells of the 
second framework with net and non-net values according to the estimated net-to-gross 
values of the first framework. 
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[0069] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: populating some or all of the cells of the 
second firamework with one or more reservoir property values. 

[0070] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, ftulher includes: populating some or all of the cells of the 
second framework with one or more porosity values. 

[0071] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: populating some or all of the cells of the 
second framework with one or more permeability values. 

[0072] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: populating some or all of the cells of the 
second framework with one or more water saturation values. 

[0073] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes: populating some or all of the cells of the 
second framework with one or more reservoir property values to generate a reservoir 
cell model; and performing a flow simulation on the reservoir cell model to generate 
one or more effective reservoir property values for the first framework. 

[0074] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes: populating some or all of the cells of the 
second framework with one or more reservoir property values to generate a reservoir 
cell model; performing a flow simulation on the reservoir cell model to generate one 
or more effective reservoir property values for the first framework; and calculating the 
variability between the effective reservoir property values for the first fi-amework. 

[0075] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, further includes: populating some or all of the cells of the 
second framework with one or more reservoir property values to generate a reservoir 
cell model; performing a flow simulation on the reservoir cell model to generate one 
or more effective reservoir property values for the first framework; calculating the 
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variability between the effective reservoir property values for the first fi-amework; and 
determining whether the rate of change in the variability between the effective 
reservoir property values remains substantially the same. 

[0076] A specific embodiment of the method identified above, or of a method 
described elsewhere herein, fiirther includes: populating some or all of the cells of the 
second firamework with one or more reservoir property values to generate a reservoir 
cell model; performing a flow simulation on the reservoir cell model to generate one 
or more effective reservoir property values for the first firamework; and populating the 
first firamework with the effective reservoir property values to generate the reservoir 
model. 

[0077] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the reservoir model is a flow simulation model. 

[0078] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the reservoir model is a geologic model. 

[0079] In at least one specific embodiment, a method of generating a reservoir 
model, includes: providing a first firamework having a plurality of cells, wherein the 
first firamework is a reservoir firamework; and providing a second firamework having a 
plurality of cells, wherein each one of the cells of the second firamework is 
substantially the same size as a sample of well data. 

[0080] In at least one specific embodiment, a method of generating a reservoir 
model, includes: providing a firamework having a plurality of cells, wherein each cell 
is the substantially same size as the well data; identifying some or all of the cells of 
the firamework as net or non-net; populating some or all of the cells of the firamework 
with one or more reservoir properties to provide a reservoir cell model; and 
performing a flow simulation on the reservoir cell model to generate one or more 
effective reservoir property values. 
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[0081] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the firamework is substantially the same size as one cell of 
a reservoir firamework. 

[0082] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the firamework is greater than the size of one cell of a 
reservoir framework. 

[0083] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, identifying some or all of the cells includes populating 
some or all of the cells of the framework with net and non-net values that correspond 
to one or more estimated net-to-gross values of a reservoir framework for the 
reservoir model. 

[0084] In a specific embodiment of the method identified above, or of a method 
described elsewhere herein, the sample of well data is the same size as a sample of 
core data. 

[0085] In a specific embodhnent of the method identified above, or of a method 
described elsewhere herein, the sample of well data is the same size as a sample of log 
data. 

Specific Embodiments In Drawings 

[0086] Specific embodiments shown in the drawings will now be described, 

[0087] Figure 1 illustrates a flow diagram of a method 100 for generating one or 
more effective reservoir properties in accordance with one embodiment of the 
invention. At step 10, a set of source data is received from the user. In one 
embodiment, such data includes a set of estimated net-to-gross values in the target 
reservour framework, which will be described in the following paragraph below. A 
detailed description of the net-to-gross is provided in the definition section of this 
application. The set of estimated net-to-gross values may be calculated using 
conventional algorithms generally known by persons of ordinary skill in the art. 
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[0088] Another set of source data that the user may specify is well data, from 
which porosity, permeability and water saturation values may be obtained. The 
significance of porosity, permeability and water saturation values will be described 
later in the following paragraphs, in particular, with reference to step 50. Well data 
includes, but not limited to, core data and log data. A more detailed description of 
well data, core data and log data is provided in the definition section of this 
supplication. 

[0089] In addition to receiving the abovementioned data, in accordance to one 
embodiment of the invention, certain modeling parameters may also be received from 
the user. Examples of those modeUng parameters include the size for each cell in a 
target reservoir framework, the size for each cell in the cell framework and various 
geologic property modeling parameters. The target reservoir framework is defined 
herein as the reservoir framework into which the effective reservoir properties are 
populated. In one embodiment, the effective reservoir properties may be populated 
into the target reservoir framework to build a target reservoir model at which flow 
simulation can be performed, i.e., a flow simulation model. In another embodiment, 
the effective reservoir properties may be populated into the target reservoir 
framework to build a target reservoir model at a geologic scale, i.e., a geologic model. 

[0090] The size for each cell in the target reservoir framework determines the size 
of the target reservoir model to be built. In other words, the size for each cell in the 
target reservoir framework determines whether the effective reservoir properties will 
be used to build a flow simulation model or a geologic model. The size for each cell 
in the target reservoir framework therefore depends upon what the user desires. 

[0091] The cell framework comprises a plurality of cells, and each cell of the cell 
framework is configured to be substantially the same size as a sample of the well data. 
In one embodiment, each cell of the cell firamework is substantially the same size as a 
sample of the core data. In another embodiment, each cell of the cell framework is 
substantially the same size as the log data sample size. The sample size of well data, 
including log data and core data, is conventional, as conunonly known by persons of 
ordinary skill in the art. 
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[0092] In addition to the modeling parameters described above, various geologic 
property modeling parameters may be provided since these parameters may affect the 
geostatistical estimation algorithm used in populating the reservoir property values, 
which will be described in the following paragr2q>hs with reference to steps 40 and 50. 
A description and/or definition for geologic modeling parameters is provided in the 
definition section above. 

[0093] At step 20, a cell fi^amework is built. The cell fi-amework may be three 
dimensional. In another embodiment, the cell framework is substantially the same 
size as one cell of the target reservoir fi:'amework. In one embodiment, the size of the 
framework is determined by the cell size of the target reservoir framework received 
from the user at step 10. As such, the cell framework is configured to be 
substantially the same size as one cell of the target reservoir framework. 

[0094] An embodiment of the cell framework is illustrated in Figure 2 as cell 
framework 200. The cell framework 200 is composed of a plurality of cells 210. In 
one embodiment, each cell 210 is substantially the same size as a sample of the well 
data received at step 10. For example, each cell 210 may be substantially the same 
size as a sample of the core data. For another example, each cell 210 is substantially 
the same size as a sample of the log data. In another embodiment, the cell firamework 
200 and each cell 210 contained therein is cubic in shape. However, the cell 
framework and the cells contained therein may be in any shape conventionally known 
by persons of ordinary skill in the art. A more detailed description of the cell 
framework and the cells contained therein is provided in the definition section of this 
application. 

[0095] Referring back to Figure 1, once the cell framework is built, a net-to-gross 
value is selected fix>m the set of estimated net-to-gross values (step 30). In one 
embodiment, the net-to-gross value may be randomly selected. The net-to-gross 
value may be in the form of a percentage. In one embodiment, the net-to-gross value 
is selected from a cumulative distribution fimction created from the set of estimated 
net-to-gross values. The cumulative distribution fimction may be derived from 
seismic data, regional maps, hand drawn maps, conceptual models, or even well- 
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based models. Alternatively, the cumulative distribution function may be created 
using conventional techniques generally known by persons of ordinary skill in the art. 

[0096] At step 40, the cells, e.g., cells 210, of the cell framework are populated 
with rock-type values that correspond to the selected net-to-gross value. In this 
manner, the cells of the cell framework are identified as rock-type 1 or rock-type 2. 
As such, one example for rock-type 1 is sand and one example for rock-type 2 is 
shale. A more detailed description of the relationship between rock-type and 
net-to-gross is provided in the definition section of this application. The cells of the 
cell framework may be populated using any conventional geostatistical estimation 
algorithm commonly known by persons of ordinary skill in the art. A more detailed 
description of geostatistical estimation algorithm is provided in the definition section 
of this application. 

[0097] At step 50, the cells within the cell framework are populated with reservoir 
property values, such as porosity and permeability, to build a reservoir cell model, 
which may be a three dimensional reservoir cell model. In one embodiment, the 
porosity and permeability values may be populated in a sequential manner. That is, 
the cells within the cell firamework may be populated with porosity values first, 
followed by permeability values. In another embodiment, the cells within the cell 
framework may also be populated with water saturation values to build the three 
dimensional reservoir framework cell model. The porosity, permeability and water 
saturation values may be obtained from the well data received from the user at step 
10, as described above. As in step 40, the porosity, permeability and water saturation 
values may be populated using any conventional geostatistical estimation algorithm 
commonly known by persons with ordinary skill in the art. Porosity, permeability and 
water saturation values are fiirther defined in the definition section of this application. 

[0098] At step 60, a flow simulation on the three dimensional reservoir cell model 
is performed to generate an effective reservoir property value, which may include an 
effective porosity value, an effective permeability value, an effective net-to-gross 
value, an effective water saturation value and an effective endpoint saturation value. 
The effective reservoir property value is configvu-ed to be populated into the cells of 
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the target reservoir framework. Effective porosity is defined as the volume- weighted 
arithmetic average of the porosity values in the target reservoir firamework cell. 
Effective permeability of a volume of a reservoir model cell is defined as a constant 
permeability value of an equivalent volume that would give the same amount of flow 
across a given area per unit time, for the same boxmdary conditions. Effective 
permeability value may include a full permeability tensor, which typically has 3 
difTerent components that are defined by 9 different values, i.e., Kx, Ky, Kz, Kxy, 
Kxz, Kyz, Kyx, Kzx and Kzy. 

[0099] When the cell framework is built at step 20, the volume of the cell 
firamework is larger than the size of one cell of the target reservoir framework. That 
is, the cell framework is buih with expanded boundary conditions, which allows a 
flow simulation to be performed without being limited by no-flow boundary 
conditions. As such, a flow simulation performed on the three dimensional reservoir 
cell model built with an expanded boundary condition provides a more realistic 
simulation pathways through the reservoir, thereby leading to more accurate effective 
reservoir properties. The size of the expanded boundary condition may be specified 
by the user at step 10 as a modeling parameter. 

[00100] At step 70, the variability between the effective reservoir property values 
is calculated. The variability may be calculated using any conventional algorithm 
known by persons with ordinary skill in the art. A more detailed description of how 
the variability between the effective reservoir property values is calculated is provided 
in the definition section of this application. 

[00101] At step 80, the rate of change in the variability between the effective 
reservoir property values is determined. If the rate of change is not near zero, then 
processing continues to step 30, at which another value from the set of estimated net- 
to-gross values is selected. On the other hand, if the rate of change is near zero, then 
processing ends. A near zero rate of change indicates that the variability between the 
effective reservoir property values has remained substantially the same. 
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[00102] In this manner, method 100 may be used to generate a set of effective 
reservoir property values, which may be populated into some or all of the cells within 
the target reservoir framework to build the target reservoir model. As such, various 
embodiments of the invention may be used as part of a method of building a reservoir 
model, as commonly known by persons with ordinary skill in the art. If the reservoir 
model is a geologic model, then the geologic model may be upscaled, e.g., using the 
technique described in commonly assigned publication WO 00/79423, published on 
December 28, 2000, and the technique described therein is incorporated herein by 
reference. Altematively, the geologic model may be upscaled prior to flow simulation 
according to conventional techniques commonly known by persons of ordinary skill 
in the art. 

[00103] Figure 3 illustrates a flow diagram of a method 300 for generating one or 
more reservoir property values in accordance with another embodiment of the 
invention. At step 3 10, a set of source data is received from the user. Such source 
data may include a set of estimated net-to-gross values in a target reservoir framework 
and well data. In addition to receiving the source data, certain modeling parameters 
may also be received from the user. Such modeling parameters may include the 
number of target reservoir framework cells to be contained inside the cell fiBmework, 
the size for each cell in the target reservoir framework, and the size for each cell in the 
cell framework. Since step 310 perforais the same function as step 10 described in 
Figure 1, the reader is directed to the paragraphs above in connection with step 10 for 
a more detailed description of step 310. 

[00104] At step 320, a three dimensional cell framework is built. The cell 
framework is configured to contain a number of target reservoir framework cells. 
This number may be provided by the user at step 310 above as one of the modeling 
parameters. In one embodiment, the number of target reservoir framework cells 
ranges from about 4 to 10. Consequently, the cell framework is much larger than the 
size of one target reservoir framework cell. As a point of distinction between step 20 
and step 320, the cell framework built at step 320 is much larger than the cell 
framework built at step 20 since the cell framework built at step 20 is substantially the 
same size as one target reservoir framework cell. Like the cell framework described 
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in step 20, the cell framework in step 320 is composed of a plurality of cells. Each 
cell is substantially the same size as a sample of the well data received at step 310. In 
one embodiment, each cell is substantially the same size as a sample of the core data. 
In another embodiment, each cell is substantially the same size as a sample of the log 
data. A more detailed description of the cell framework and the cells contained 
therein is provided above in connection with step 10. 

[00105] Processing then continues to steps 330 through 350. However, steps 330 
through 350 perform the same ftmctions as steps 30 through 50. Accordingly, the 
reader is directed to the paragraphs above in connection with steps 30 through 50 for a 
detailed description of steps 330 through 350. 

[00106] At step 360, a cell sample is randomly extracted from the three 
dimensional reservoir cell model. In one embodiment, the cell sample is substantially 
the same size as one target reservoir framework cell. In another embodiment, the cell 
size is larger than one target reservoir framework cell to allow for expanded boundary 
conditions (see paragraph 98 above). 

[00107] At step 370, a flow simulation on the cell sample is performed to generate 
an effective reservoir property value, which may include an effective porosity value, 
an effective permeability value, an effective net-to-gross value, an effective water 
saturation value and an effective endpoint saturation value. A more detailed 
description of the effective reservoir property value, an effective porosity value, an 
effective permeability value, an effective net-to-gross value, an effective water 
saturation value and an effective endpoint saturation value is provided above with 
reference to step 60. 

[00108] Processing then continues to step 380, which performs the same function 
as step 70. Accordingly, the reader is directed to the paragraphs above in connection 
with step 70 for a detailed description of step 380. 

[00109] At step 390, the rate of change in the variability between the effective 
reservoir property values is determined. If the rate of change is not zero, then 
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processing continues to step 395, which will be described in the following paragraph 
below. On the other hand, if the rate of change is zero, then processing ends. 

[00110] At step 395, a determination is made as to whether a sample of the three 
dimensional reservoir framework cell model has been extracted according to the 
number of target reservoir framework cells specified by the user at step 310. If the 
answer is in the negative, then processing returns to step 360, at which another sample 
of the target reservoir framework cell is extracted from the three dimensional 
reservoir fiBmework cell model. If the answer is in the affirmative, then processing 
returns to step 330, at which another net-to-gross value is selected from the set of 
estimated net-to-gross values. In this manner, a number (as specified by the user at 
step 330) of samples of the three dimensional reservoir framework cell model is 
sampled for each net-to-gross value selected at step 330. 

[001 111 Embodiments of the invention have many advantages. For example, 
embodiments of the invention eliminate the need to building a large, fiiU-field 
reservoir model and upscaling the reservoir model. For another example, 
embodiments of the invention use a statistical sampling procedure to minimize the 
number of fine-scale simulations required and provide libraries of effective reservoir 
properties for specific geologic features, which can be maintained and used for other 
reservoirs. Further, embodiments of the invention provide a better statistical 
treatment of core data accoxmting for the volume differences between core plugs and 
geologic model cells. 

[00112] Figure 4 illustrates a computer network 400, into which embodiments of 
the invention may be implemented. The computer network 400 includes a system 
computer 430, which may be implemented as any conventional personal computer or 
workstation, such as a UNIX-based workstation. The system computer 430 is in 
communication with disk storage devices 429, 43 1, and 433, which may be external 
hard disk storage devices. It is contemplated that disk storage devices 429, 431, and 
433 are conventional hard disk drives, and as such, will be implemented by way of a 
local area network or by remote access. Of course, while disk storage devices 429, 
431, and 433 are illustrated as separate devices, a single disk storage device may be 
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used to store any and all of the program instructions, measurement data, and results as 
desired. 

[001 131 ^ embodiment, the input data are stored in disk storage device 43 1 . 
The system computer 430 may retrieve the appropriate data from the disk storage 
device 43 1 to perform the reservoir model generation according to program 
instructions that correspond to the methods described herein. The program 
instructions may be written in a computer programming language, such as C-+H-, Java 
and the like. The program instructions may be stored in a computer-readable 
memory, such as program disk storage device 433. Of course, the memory medium 
storing the program instructions may be of any conventional type used for the storage 
of computer programs, including hard disk drives, floppy disks, CD-ROMs and other 
optical media, magnetic tape, and the like. 

[00114] According to a preferred embodiment, the system computer 430 presents 
output primarily onto graphics display 427, or alternatively via printer 428. The 
system computer 230 may store the results of the methods described above on disk 
storage 429, for later use and further analysis. The keyboard 426 and the pointing 
device (e.g., a mouse, trackball, or the like) 225 may be provided with the system 
computer 430 to enable interactive operation. 

[001151 The system computer 430 may be located at a data center remote from the 
reservoir. While Figure 4 illustrates the disk storage 431 as directly connected to the 
system computer 430, it is also contemplated that the disk storage device 431 may be 
accessible through a local area network or by remote access. Furthermore, while disk 
storage devices 429, 431 are illustrated as separate devices for storing input data and 
analysis results, the disk storage devices 429, 431 may be implemented within a 
single disk drive (either together with or separately from program disk storage 
device 433), or in any other conventional manner as will be fully understood by one 
of skill in the art having reference to this specification. 



